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WHY MASSIVE BLACK HOLES ARE SMALL IN DISK GALAXIES ? 
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ABSTRACT 

A potential mechanism is proposed to account for the fact that supermassive black holes (SMBHs) 
in disk galaxies appear to be smaller than those in elliptical galaxies in the same luminosity range. We 
consider the formation of SMBHs by radiation drag (Poynting- Robertson effect), which extracts angular 
momentum from interstellar matter and thereby drives the mass accretion onto a galactic center. In 
particular, we quantitatively scrutinize the efficiency of radiation drag in a galaxy composed of bulge 
and disk, and elucidate the relation between the final mass of SMBH and the bulge-to-disk ratio of the 
galaxy. As a result, it is found that the BH-to-galaxy mass ratio, AfEnZ-^^gaiaxy, decreases with a smaller 
bulge-to-disk ratio, due to the effects of geometrical dilution and opacity, and is reduced maximally by 
two orders of magnitude, resulting in MBH/-Wgaiaxy ~ 10~^. In contrast, if only the bulge components 
in galaxies are focused, the BH-to-bulge mass ratio becomes MeH/Afbuigc ~ 10~^, which is similar to 
that found in elliptical galaxies. Thus, it turns out that the mass of SMBH primarily correlates with a 
bulge, not with a disk, consistently with observational data. 

Subject headings: black hole physics — galaxies:active — galaxies:disk — galaxiesmuclei — 
galaxies:starburst 



1. INTRODUCTION 

Supermassive black holes (SMBHs) are now thought to 
reside in virtually all galactic bulges. The recent compi- 
lation of the kinematical data of galactic centers in both 
active and inactive galaxies has shown that the mass of 
SMBH does correlate with the mass of a galactic bulge; 
the mass ratio of BH to bulge is « 0.001 as a median 
value (Kormendy & Richstone 1995; Richstone et al. 1998; 
Magorrian et al. 1998; Gebhardt et al. 2000; Ferrarese & 
Merritt 2000; Merritt & Ferrarese 2001; McLure & Dunlop 

2001, 2002; Wandel 2002). On the other hand, it appears 
that the SMBH does not correlate with a disk component 
(Kormendy & Gebhardt 2001). In disk galaxies, the mass 
ratio is significantly smaller than 0.001 if disk stars are 
included (Salucci et al.2000; Sarzi et al. 2001, Ferrarese 

2002, Baes et al. 2003). These observations imply that the 
formation of the SMBH is related to the bulge components 
of galaxies. 

Recently, as a possible mechanism for the formation of 
SMBHs, the radiation-hydrodynamic model is proposed, 
where the radiation drag (Foynting-Robertson effect) due 
to stellar radiation in a galactic bulge drives the mass 
accretion to grow a massive BH (Umemura 2001). This 
is a parallel model of the formation of massive BHs by 
the Compton drag in an early universe (Umemura, Loeb, 
& Turner 1997). In an optically thick regime, the efh- 
ciency of radiation drag is saturated due to the conserva- 
tion of the photon number (Tsuribe & Umemura 1997). 
The timescale of radiation drag-induced mass accretion is 
idrag = 8.6 X lO7i?2p^(i,/lOi2Lo)-i(^/^0)"Vr, where 
^kpc is the galaxy size in units of kpc, and Z is the 
metallicity of gas. The mass accretion rate is given by 
M = —AlgdlnJ/dt ~ L^./c^, where J, L*, and Mg are the 
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total angular momentum of gaseous component, the total 
luminosity of the system, and the total mass of gas, respec- 
tively. Thus, the final BH mass is estimated by Mbh — 
J L^,{t)/c^dt. Hence, the final BH mass is proportional to 
the total radiation energy from stars, and the theoretical 
upper limit of BH-to-bulge mass ratio is determined by the 
energy conversion efficiency of nuclear fusion from hydro- 
gen to helium, i.e., 0.007 (Umemura 2001). In practice, 
the inhomogeneity of interstellar medium helps the radia- 
tion drag to sustain high efficiency (Kawakatu & Umemura 
2002). By incorporating the realistic chemical evolution, 
the radiation drag model predicts MBH/A/buigc — 0.001 
(Kawakatu, Umemura, & Mori 2003). However, the radia- 
tion drag efficiency could be strongly subject to the effect 
of geometry (Umemura et al. 1997,1998; Ohsuga et al. 
1999a). If the system is spherical, the emitted photons are 
effectively consumed within the system, whereas a large 
fraction of photons can escape from disk- like systems, and 
thus the drag efficiency is likely to be reduced considerably. 
It is suggested that such a geometrical effect may be the 
reason why the BH is smaller in disk galaxies (Umemura 
2001; Umemura 2003), although the details have not been 
clear quantitatively. 

In this Letter, we scrutinize the efficiency of radiation 
drag in a system composed of bulge and disk. For the 
purpose, we accurately solve the 3D radiative transfer in 
an inhomogeneous interstellar medium. The main goal is 
to elucidate the relationship between the final BH mass 
formed by the radiation drag and the bulge-to-disk ratio 
of host galaxy. 

2. MODEL 

We consider a system composed of bulge and disk, which 
consist of dark matter, stars, and dusty interstellar mat- 
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ter (ISM). The total baryonic mass of the system, Mb, 
is set to be 10^^ Mq. But, it is noted that the results 
are scaled with Mb. The dark matter is distributed fol- 
lowing the profile by Navarro, Frenk, & White (1997), 
p(r) cx (r/rs)^^(l + r/rs)~^, where is the character- 
istic radius, which is Vg = 50kpc here. The mass of the 
dark matter component, Mdm, is equal to Mb within the 
galaxy size of Vg =10kpc. 

We incorporate the evolution of galaxy from the for- 
mation stage, where a galactic systcini starts with a pro- 
togalactic cloud purely composed of barynic gas. By in- 
voking the star formation in Schmidt's law, stellar com- 
ponents are generated in this system. To include the re- 
cycling of ISM, we employ an evolutionary spectral syn- 
thesis code, PEGASE (Fioc & Rocea-Volmerange 1997), 
with assuming a Salpeter-type initial mass function (IMF) 
as (/> = dn/dlogm* = (?!>o(m,/MQ)~-^-^^ for the mass range 
of [O.lAfQ, 6OM0], where is the stellar mass. The 
star formation rate (SFR) per unit mass at time t, C'(t), 
is assumed to be proportional to the gas mass fraction, 
C{t) = kfg with /g = Mg(i)/Mb, where /c is a rate co- 
efficient which is 8.6Gyr^^ for the bulge component and 
0.86Gyr~^ for the disk component, respectively. 

The formed stars in the bulge are distributed following 
Hernquist's profile (1990), which mimics de Vaucouleur's 
profile: p{r) cx (r/rcorc)~Hl + ^Acoro)"^ , where rcoro is 
the core radius. The bulge radius, rbuige; is given by a 
mass-to-size relation as rbuige /rg = (Mbulge/Mgalaxy)^/^ 
(Kormendy 1985; Chiosi & Carraro 2002), where Mbuigc 
is the bulge mass and Mgaiaxy is the final stellar mass of 
the galaxy, and rcore is set to Tcore = O.lrbuige- In the 
present analysis, the bulge mass fraction (Mbuige/^fgaiaxy) 
is altered as a key parameter. In addition, angular mo- 
mentum is smeared, assuming rigid rotation, according to 
the spin parameter of A = 0.05 (Barnes & Efstathiou 1987; 
Heavens & Peacock 1988). 

As for the disk component, formed stars are distributed 
exponentially in the radial direction, with assuming the 
vertical scale-height which is changed from O.Olrg to O.lrg. 
The disk radius, rdisk) is given by a mass-to-size relation 
of rdisk/rg = (Mdisk/Afgaiaxy)^/^, where Mdisk is the disk 
mass. The disk component is assumed to be in rotation 
balance. 

As for ISM, we consider clumpy matter, since the ISM 
is observed to be highly inhomogeneous in active star- 
forming galaxies (Sanders et al. 1988; Gordon, Calzetti, 
& Witt 1997). Here, Nci= 10'*) identical clouds are dis- 
tributed randomly in a system of the bulge and disk. The 
number of clouds is minimal enough to resolve the disk 
seal-height of O.Olrg. (It is noted that simulations with a 
three times larger number of clouds did not lead to any 
fundamental difference.) The distribution and motion of 
the ISM is the same as those of stellar components. We 
assume that the cloud covering factor is unity according 
to the analysis by Kawakatu & Umemura (2002). The in- 
ternal density in a cloud is assumed to be uniform, and 
the opacity is determined by a dust-to-gas ratio which is 
calculated by PEGASE. The optical depth of a cloud and 
therefore the overall optical depth of galaxy, tt, depend 
on the cloud size rc, where Tc = lOOpc is assumed as a 
fiducial case. But we have confirmed that no essential dif- 
ference in the final BH mass is found by changing rc so as 



to enhance tt by an order. 
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Fig. 1. — The BH-to-galaxy mass ratio (Mbh /JWgaiaxy ) 
and the BH-to-bulge mass ratio (MBn/A^buige) are shown against 
the bulge fraction, /bulge = Afbuige /-^galaxy The thick hnes 
show MenZ-'l^galaxy with changing the disk scale- height from h = 
O.Olrdisk up to /i = O.lrdisk and the thin lines do MBu/Afbulge- The 
top abscissa is the indication for the Hubble type, which is deter- 
mined by the mean bulge-to-disk ratio for disk galaxies (Kent 1985) . 

The radiation drag, which drives the mass accretion, 
originates in the relativistie effect in absorption and ree- 
mission of radiation. This effect is involved in relativistie 
radiation hydrodynamics (RHD) (Umemura, Fukue, & Mi- 
neshige 1997, 1998; Fukue, Umemura, & Mineshige 1997). 
The angular momentum transfer in RHD is given by the 
azimuthal equation of motion in cylindrical coordinates, 

ld{n^^X^^^_ (1) 
r at c 

where x is the mass extinction coefficient, E is the radi- 
ation energy density, F'^ is the radiation flux, and P'f"f' 
is the radiation stress tensor. All radiative quantities are 
determined by radiation from stars diluted by dusty ISM. 
Here, we solve radiative transfer including the dust extinc- 
tion by ISM and obtain these radiative quantities. Then, 
we evaluate the total angular momentum loss rate (see 
Kawakatu & Umemura 2002 for the details of method). 
If an optically thick cloud is irradiated by the radiation, 
an optically-thin surface layc;r is stripped by the radia- 
tion drag (Tsuribe & Umemura 1997), and loses angular 
momentum. Since the bulge luminosity is lower than the 
Eddington luminosity for dusty gas, LEdd = 47rcGMdyn/x 
with the dynamical mass M^yn within rg, during the whole 
history of galaxy, the stripped gas is not blown away by 
the radiation pressure, but confined in the galaxy. Hence, 
the gas moves in a gravitational potential weakened by 
radiation pressure and loses angular momentum by radi- 
ation drag. The stripped gas losing angular momentum 
falls onto the galactic center to form a massive object, 
which is likely to undergo the inside-out viscous collapse 
(Mineshige, Tsuribe, & Umemura 1998; Tsuribe 1999) and 
eventually evolve into a SMBH by the general relativistie 
instability (Saijo et al. 2002; Shibata & Shapiro 2002). 
Using equation (1), we can calculate the mass accretion 
rate onto the central object by Mg/Mg = —J/ J, where 
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J is the total angular momentum of ISM. Then, the BH 
mass (Mbh) is assessed as 

Mbh = / Mgdt = - M^-dt, (2) 
Jo Jo J 

where to is the Hubble time. 
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Fig. 2. — Comparison between the theoretical prediction and the 
observational data on BH mass for disk galaxies. The upper panel 
shows the BH-to-galaxy mass ratio (MBn/^galaxy) ^nd the lower 
panel shows the BH-to-bulge mass ratio (MBHZ-'W^bulgc) against the 
bulge fraction, /bulge = A^buige/A^gaiaxy The observational data 
are based upon Table 1 and plotted by symbols: Crosses - disk 
galaxies which do not possess AGNs, open circles - Seyfert 1 galax- 
ies (Syls), filled triangles - narrow line Seyfert 1 galaxies (NLSyls), 
and filled circles - Seyfert 2 galaxies (Sy2s). Seyfert galaxies ac- 
companied by starburst activities are specified like Syl / starburst or 
Sy2/starburst. The hatched area is the prediction shown in Fig. 1. 
Two dotted lines in each panel are the upper and lower bound of 
the theoretical prediction by including uncertainties other than the 
geometrical effect (see the end of §3). 

3. RESULTS 

By changing the bulge mass fraction, /bulge = 
.^buige/.Wgaiaxy, wc investigate the relation between the 
final BH mass and the bulgc-to-disk ratio of galaxy. The 
growth of BH depends upon the optical depth of galaxy, 
tt, but no essential difference in the final BH mass is found 
by enhancing tt by an order, because the mass accretion 
rate is saturated when the system becomes optically-thick 
(Umcmura 2001; Kawakatu & Umcmura 2002). In Figure 
1, the resultant BH-to-galaxy mass ratio, MeH/Afgaiaxy, 
is shown against the bulge mass fraction (/bulge) by thick 
lines for three values of disk scale-height. We can see a 
clear trend that MeH/Afgaiaxy decreases with decreasing 
/bulge- This tendency comes from three effects: First, a 
larger number of photons escape from the disk surface, 
since the surface-to- volume ratio is larger for the disk com- 



ponent compared to the bulge. This effect is called the 176- 
ometrical dilution here. Secondly, the radiation from disk 
stars is heavily diminished across the disk, because the 
edge-on optical depth becomes larger quickly. Thirdly, 
the velocity difference between a star and an absorbing 
cloud becomes closer to zero in an optically-thick disk, so 
that the radiation drag cannot work effectively. For in- 
stance, if /bulge = 0.03 which corresponds to Sd galaxies, 
M^Yi/Mgaia^y IS rcduccd by a factor of « 20 to ~ 100, 
compared to MeH/A/gaiaxy ~ 0.001 for elliptical galaxies 
(/bulge = !)• Furthermore, it is predicted that a bulgeless 
(pure disk) system can also possess a massive black hole 
with the mass ratio of MBH/-/^gaiaxy — 0.001 x (/i/2rdisk), 
where h is the scale height of disk. This ranges from 
6 X 10^^ to 5 X 10~^, corresponding to h = O.Olrdisk and 
h = O.lrdisk, respectively. 

On the other hand, the bulge component itself is likely 
to be impervious to the above three geometrical effects. 
In Figure 1, by focusing on only the bulge component, 
the BH-to-bulge mass ratio, MBH/.Wbuige) is also shown. 
As seen in Figure 1, MBn/Afbuige is almost constant as ex- 
pected. Hence, it is concluded that the BH-to-galaxy mass 
ratio in disk galaxies is substantially determined by bulge 
components, not by disk components: MBH/-Mgaiaxy — 
0.001, where the constant value of 0.001 is basically deter- 
mined by the energy conversion efficiency of nuclear fusion 
from hydrogen to helium, i.e., 0.007 (Umemura 2001). The 
BH accretion can trigger an active galactic nucleus (AGN) 
as well, and then the further mass accretion is induced 
via radiation drag by the AGN. The effect can enhance 
the final BH mass maximally by a factor of 1.7 (Umemura 

2001) . Also, a different level of ISM dumpiness can reduce 
the final BH mass by a factor of 2 (Kawakatu & Umemura 

2002) . As for the effect of IMF, if the slope and the mass 
range of IMF are changed to satisfy the spectrophotomet- 
ric properties, then the final BH mass is altered by a factor 
of ±40% . Similar uncertainties attach also to the final BH 
mass in disk galaxies. 

4. COMPARISON WITH OBSERVATIONS 

The present results are compared with observational 
data for SMBHs in disk galaxies in Figure 2. The data we 
adopted are summarized in Table 1. In the upper panel 
of Figure 2, MBn/Mgaiaxy is shown against /bulge- The 
hatched area is the prediction shown in Figure 1. Two 
dotted lines are respectively the upper and lower bound 
of the prediction by including other theoretical uncertain- 
ties mentioned above. The observational data are cate- 
gorized into four types: disk galaxies which do not pos- 
sess AGNs, Seyfert 1 galaxies (Syls), narrow line Seyfert 
1 galaxies (NLSyls), and Seyfert 2 galaxies (Sy2s). Also, 
the objects accompanied with starburst events are speci- 
fied. The observational data in Figure 2 shows the decrease 
of MBH/-M^gaiaxy by morc than an order of magnitude from 
/bulge = 1 to /bulge — 0.03. This trend is broadly consis- 
tent with the theoretical prediction. But, it is worth noting 
that several objects possess relatively small BHs compared 
to the prediction, and they are mostly Seyfert galaxies 
with starburst events (NGC 7469, NGC 1068, Circinus) or 
NLSyls (Mrk 590, NGC4051). In recent years, an evolu- 
tionary model for AGNs beyond the AGN unified theory 
has been considered from both observational and theo- 
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rctical points of view (Maiolino et al. 1995; Radovich, 
Rafanelli, & Barbon 1998; Ohsuga & Umemura 1999, 
2001). In this model, AGNs with starbursts are in earlier 
evolutionary phase, and therefore the BH mass is likely to 
be smaller, compared to a post-starburst phase. Further- 
more, NLSyls are thought to be in a rapidly growing phase 
of an immature BH (Mathur et al. 2001). NGC1068 could 
be also in such a phase (Kawaguchi 2003). This point may 
be intriguing to analyze further with larger sample data 
in the future. 

In the lower panel of Figure 2, MsH/Afbuige is shown 
against /bulge- When only the bulge component is focused, 
the observational data for MBH/-/^buigo lies at a level of 
w 0.001, and they roughly agree with the prediction. But, 
again the ratios for Seyferts with starbursts and NLSyls 



fall appreciably below 0.001. 

The present model also predicts a SMBH even in a pure 
disk system without bulge, and the mass of BH depends 
upon the disk scale-height. NGC4395 is an archetype of 
bulgclcss galaxy in which a massive BH resides (Filipcnko 
& Ho 2003). The BH fraction in this galaxy is estimated to 
be MBH/-^gaiaxy — 1-5 X 10~^. In the present model, this 
level can be realized by a disk scale height of /i « 0.04rdisk- 

We thank T.Nakamoto, T.Saitou and A.Yonehara for 
helpful discussions. We also thank the anonymous ref- 
eree for precious comments. Numerical simulations were 
performed with facilities at the Center of Computational 
Physics, University of Tsukuba. This work was supported 
in part by the Grant-in- Aid of the JSPS, 11640225. 
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Tabic 1. SMBHs in Disk Galaxies 



Object 


Type 




xy(M0) 




Mbh{M 






AGn 


Starburst 


Reference 


Galaxy 


Sbc 


4.6 X 


lO"' 


5.2 X 10'' 


3.7(±0.4) 


X 


lO*" 






1 


M31 


Sb 


1.4 X 


IQii 


1.9 X 10^" 


7(+13 - 4) 


X 


10^ 






2 


M81 


Sb 


9.7 X 


10^° 


1.8 X 10^" 


6.8(+0.7 - 1.3) 


X 


10^ 






3 


NGC1023 


SBO 


1.4 X 


10^1 


2.3 X 10^° 


4.4(+0.4 - 0.5) 


X 


lO'' 






4 


NGC3245 


SO 


7.5 X 


lOi" 


3.9 X 10^" 


2.1(±0.5) 


X 


10** 






5 


NGC7457 


SO 


2.1 X 


lOi" 


9.3 X lO** 


3.5(+l.l - 1.4) 


X 


10'' 






6 


Fairall 9 


S 


2.2 X 


10" 


2.9 X 10^" 


8.7(+2.6 - 4.5) 


X 


10^ 


Syl 




7 


Mrk509 


s 


2.9 X 


10" 


5.7 X 10^" 


9.5(+l.l - 1.1) 


X 


10^ 


Syl 




7 


NGC3516 


SBO 


8.0 X 


10^° 


4.9 X 10" 


1.7(±0.3) 


X 


10^ 


Syl 




8 


NGC3783 


SBa 


1.4 X 


IQio 


4.9 X 10^ 


1.1(+1.1 - 1.0) 


X 


10^ 


Syl 




7 


NGC4151 


SA 


5.7 X 


10» 


2.0 X 10^ 


1.2(+0.8 - 0.7) 


X 


10^ 


Syl 




7 


NGC3227 


SAa 


5.6 X 


10« 


2.9 X 10^ 


3.6(±1.4) 


X 


10^ 


Syl 




8 


NGC4593 


SBb 


3.3 X 


10^" 


1.6 X 10" 


6.6(±5.2) 


X 


10** 


Syl 




8 


NGC5548 


Sa 


4.4 X 


lO^o 


2.0 X 10" 


8.0(+1.5 - 1.0) 


X 


10^ 


Syl 




7 


NGC7469 


SA 


8.9 X 


lO^" 


5.5 X lO" 


7.6(+7.5 - 7.6) 


X 


10** 


Syl 


Nuclear 


7 


3C120 


SO 


2.2 X 


10" 


5.2 X 10^" 


3.0( + 2.0 - 1.5) 


X 


10^ 


NLSyl 




7 


Mrk590 


Sa 


2.3 X 


10" 


1.1 X lO" 


1.4(+0.3 - 0.3) 


X 


10^ 


NLSyl 




7 


NGC4051 


SB 


5.6 X 


10^" 


1.1 X lO"' 


1.4(+1.5 - 0.9) 


X 




NLSyl 




7 


Circinus 


Sb-Sd 


1.9 X 


10" 


9.0 X lO" 


1.7(±0.3) 


X 


lO" 


Sy2 


Nuclear 


9 


NGC1068 


Sb 


7.2 X 


10" 


3.6 X lO" 


1.7(+1.3 - 0.7) 


X 


10" 


Sy2 


Circumnuclear 


10 


NGC4258 


Sbc 


1.3 X 


10" 


6.2 X 10^ 


4.0(±0.1) 


X 


10" 


Sy2 




11 


NGC4945 


Scd 


2.4 X 


10^° 


7.2 X 10^ 


1.4(+0.7-0.5) 


X 


10" 


Sy2 


Nuclear 


12 
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